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1 Introduction 



This paper has two motivations. On the one hand, we are interested in physical models of 
particles on a number of half-lines joined at a central point. Such systems, for free particles, 
have been treated in, for example, [I1I2]- Here we would like to consider interacting models, to 
establish that integrable examples of such models exist, and to find their symmetries. We shall 
work out explicitly two examples: the Gaudin model [3] and the Calogero model [3]. Both have 
numerous applications in physics and in mathematics. For example, the reduced BCS model for 
conventionnal superconductivity can be diagonalized in an algebraic way [5j using the Gaudin 
model. Other, more recent, applications of the Gaudin model in quantum many-body physics 
can be found for example in the reviews OEj. Besides being of intrinsic interest due to its exact 
solvability, the Calogero model plays a role in the study of two dimensional Yang-Mills theory 
[8], the quantum Hall effect [9] and fractional statistics [10] . 

Our second motivation is algebraic. The notation of the St Peterburg school [TT] is a powerful 
tool when working with the Yangian [12] of gl^y and its subalgebras: the reflection algebras p^ [T^ 
and twisted Yangians [15]. These are quantum algebras, but the construction has a classical 
limit in which the quantum i?-matrix and Yang-Baxter Equation are replaced by their classical 
counterparts (see for example [E]). The classical limit of the Yangian is the half loop algebra, 
and the limits of the reflection algebras and twisted Yangian are subalgebras of this half loop 
algebra defined by automorphisms of order 2. But there also exist, at least in the classical case, 
other subalgebras of the half loop algebra, defined by automorphisms of higher finite order. We 
wish to study these subalgebras using classical r-matrix techniques. 

It is well-known that the half-loop algebras associated to Lie algebras are crucial in the study 
of Gaudin models and Calogero models. These algebras provide, in the former case, a systematic 
way to construct the model (see e.g pTfj) and, in the latter case, the symmetry algebras of the 
system [181 ttll 120]. In both cases, they allow one to prove the integrability of the model . We 
shall find similar connections in the cases studied in this paper. Indeed, we shall see below that 
the order n subalgebras of the half loop algebra appears naturally in the description of models 
on n half-lines. 

This paper is structured as follows. We begin with a brief review of the half loop algebra of 
gl^ and its subalgebras associated to automorphisms of order n. We make use of the notation of 
the St Petersburg school to find Abelian subalgebras. In the subsequent sections these algebraic 
results are shown to provide new quantum integrable models and demonstrate their symmetries: 
section [3] discusses "twisted" Gaudin magnets, and section U] introduces Calogero-type models 
on n half-lines joined a central point. We end with some conclusions and a short discussion of 
classical counterparts of these results. 
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2 Half loop algebra and subalgebras 



2.1 St Petersburg notation and half loop algebra 

The half loop algebra TIn based on is the complex associative unital algebra with the following 
set of generators | 1 ^ j ^ A^, a € Z^oj) subject to the defining relations 

v^ij ^^ki \- hi - "il *kj y^-^) 

for a, (3 ^ and 1 ^ i,j,k,l ^ N. It is isomorphic to the algebra fltTvl-^] of polynomials in an 
indeterminate z with coefficients in gt^y, with the generators identified as follows: 

4f = e.,^", (2.2) 

where Cij are the generators of glj^, satisfying the commutation relations 

[eij, en] = Sjkeu - duCkj. (2.3) 

It will simplify our computations to introduce the notation of the St Petersburg school: let 
Eij be the N x N matrix with a 1 in the ijth slot and zeros elsewhere. These are the generators 
of q[]\i in the fundamental representation. Let us now gather the generators of TCjy in the matrix 

N ^{a) N 

T{u) = E ^ E ^ = E ® ^d^^ = E ^ ' (2.4) 

i,j=l a^O *)i=l o^O 

where T^") = .^^ (a ^ 0) and n is a formal parameter called the spectral parameter. 

Note the flip of the indices between Eij and tji, which will prove convenient later. The algebraic 
object T{u) is an element of MatATxAf ^ 'Hn[[u^^]], and as usual we refer to MatArxAr as the 
auxiliary space and 7i]\f as the algebraic space. In what follows we shall require several copies 
of both spaces. We use letter a,b, . . . from the start of the alphabet to refer to copies of the 
auxiliary space and numerals 1,2, .. . for copies of the algebraic space. Let us introduce also 

rabiu) = ^ (2.5) 
u 

where Pat = Yli^j=i ^ij ® Eji is the permutation operator between two auxiliary spaces: the 
letters a and b stand respectively for the first and the second spaces. By definition, it satisfies 
Pab V ^ w = w CiS) V (u,wG C^). The matrix rab{u), usually called the classical R- matrix (see for 
example [IS]), satisfies the classical Yang-Baxter equation 

[rab{Ua - Ub),raciUa - life)] + [VabiUa - Ub),rbc{ub - Uc)] + [rac(Ua " Uc),rbc{ub - Uc)] = (2.6) 

and allows us to encode the half loop algebra defining relations ()2.ip in the simple equation 

[Taiu),Tbiv)] = [Ta{u)+Tb{v),rab{u-v)] . (2.7) 
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This form of commutation relations can be obtained easily by taking the classical limit of the 
presentation of the Yangian of |12j introduced by L.D. Faddeev, N.Yu. Reshetikhin and 
L.A. Takhtajan of St Petersburg [11]. By taking the trace in the space a in (|2.7p . it is straight- 
forward to show that the coefficients of the series t{u) = tr aTa{u) are central. The quotient of 
the algebra TCn by the relation t{u) = is isomorphic to the polynomial algebra 5lAr[z]. 

The identification (j2.2p between the generators of TIn and 0ljv[2:] now reads 



Taiu) = , (2.8) 

u — z 

where 

N 

P„i = ® e^-i (2-9) 

and is to be understood as the formal series Yla>o ^a+i ■ Note the similarity between relations 
(j2.5|) and (|2.8|) : the only differences are that the second auxiliary space (denoted b) in (|2.5|) is 
replaced by an algebraic space (denoted 1) and that the spectral parameter is shifted by z. In 
fact, there exists a more general solution of the relations (j2.7|) in the L-fold tensor product of 

Ta{n) = j2^- (2.10) 

From now on, we work in the enveloping algebra U{7ij\f) in which, for example, the product 
Ta{uY makes sense. 



2.2 The inner-twisted algebras 



Let a be an inner automorphism of gt^ of order n. One way to define a is by its action on 
matrices X G MatATxTv in the fundamental representation: 

a-.X^G'^XG (2.11) 

where G G MatATxAr satisfies G" = 1; the action of a on the abstract algebra gl^y is then given 
by a : Sij i— > Gjq{G~^)piepq, or, in the notation of the previous section, 

a:Fai^GaraiG-\ (2.12) 



The eigenvalues of a are the n-th roots of unity r'^ := exp 2Trik/n, and for each A; G Z„ := ^ 
the map 

Pk = -y] r-^^a^ (2.13) 
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is the projector onto the r'^-eigenspace: 

aPk = T^Pk, PkPj = SjkPk. (2.14) 

Since 1 = Pq + A + • • • + Pn-i, Q^n decomposes into the direct sum of eigenspaces of a. This 
decomposition respects the Lie bracket, in the sense that if ae = r^e and af = f then 

a[eJ] = [ae,af]=T^^\eJl (2.15) 

and is said to be a Z„-gradation of gt^v. 

By a change of basis we can take 

G = diag(l_J,Tv_r, . . . , f-\...,T^-^ ) (2.16) 

where A'^o + A^i + • • • + A''„_i = N . Note that the +l-eigenspace of a is the Lie subalgebra 
Let us define 

Q^zf = {Aiz) G gl^lz] I aAiz) = Airz)} ; (2.17) 

that is, B[Ar[2;]'^ is the subalgebra of 0lAf[-2^] which each element of degree k is also in the r^- 
eigenspace of a. There is a surjective projection map gtTvI-^] 9^n[^T j defined by ez^ i— > Pkez^ . 
In view of (|2.12|) . this sends 

T{u) ^ T^G^T{uT^)G-^ =: B{u), (2.18) 

which defines the formal series B{u) whose expansion 

B(n) = i5W + ^B«+..., (2.19) 
contains by construction a complete set of generators of 0(^[z]'^. 

Lemma 2.1 B{u) obeys 



[Ba{u),Bb{v)]= J2 

keZn 



f^—k p /^k 

r^B,{u) + B,{v),^ 



U — T^V 



(2.20) 



and has the property that for all k £ "Zn 

B{u) = t^G''B{ut^)G-^. (2.21) 

Proof. These first of these is true by virtue of (j2.7p while the second follows immediately from 
the definition (f2J8]l . ■ 

The coefficients in the expansion of b{u) = trB{u) are central in ^[^[z]"", as may be seen 
by taking the trace in space a or 6 in (|2.20p . But there exist also other abelian subalgebras in 
U{glj^[z]'^), as follows. 
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Proposition 2.2 The coefficients in the expansion ofb'{u) = trB{u)'^ are mutually commuting, 
or equivalently 

[b'{u),b'{v)] = (2.22) 
for all values of u and v. Moreover, they commute with the generators of qIj\i[z]'^ of degree zero: 

[5(0),6'(n)] =0 . (2.23) 

The algebraic elements in B^^^ generate glj^f^ © gl^^ © • • • © d'^N^-i ■ 

Proof. The details of the proof are given in appendix [XI ■ 

In particular, we recover (for n = 1) the fact that trT(M)^ commute and (for n = 2) the 
results of Hikami [19] concerning the classical limit of the reflection algebra. 

In sections [3] and [U we will apply this purely algebraic result to find new integrable models. 
2.3 Outer Automorphisms 

In the previous section, we focused on inner automorphisms. Now, we show how to modify 
the construction to study outer automorphisms. Modulo inner automorphisms, the only outer 
automorphism of gljy is generalized transposition, which has order 2. 

Let K he a real invertible N x N matrix satisfying = rjK with r/ = ±1 (for 77 = — 1, 
must be even), and define an outer automorphism T by eij ^ Kjp{K~^)giepq, or equivalently 

T : P,i ^ fJ{ = Ka^llK-^ =: Qai , (2.24) 

where ta is matrix transposition in the space a. The eigenvalues of T are ±1 and, as before, the 
decomposition of gl^v into the direct sum of eigenspaces of T defines a Z2-gradation. 

One may introduce the N x N matrices 

g+ = diag(l,...,l,-l,...,-l) and C?" = diag(l, . . . , 1) © ( ^ ) , (2.25) 




where p + q = N and the second case is valid only for N even. A well-known result in linear 
algebra is then that K is congruent over the reals to , i.e. lAKlA^ = for some real matrix 
U. From this one sees that the +l-eigenspace of T is the Lie subalgebra 5o{p, q) for rj = +1 and 
5p{N) for 7? = -1. 

Once more we may now define 

g[^[z]^ = {A{z) G Ql^[z] I TA{z) = Ai-z)} , (2.26) 
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that is, the subalgebra of Qln[z] in which each element of degree k is also in the (— l)*''-eigenspace 
of T. The projection map BlAr[2;] qIj^[z]'^ is ez'' ^ ^(1 + {—l)''T)ez'' , and, given (I2.24|) . this 
sends 

T(n) ^ r(n) + r(-u)^ =: S{u), (2.27) 
which defines the formal series S{u), whose expansion in inverse powers of u 

S{u) = -S^^^ + ^^(1) + . . . , (2.28) 

contains a complete set of generators of 0[jv[-z]'^- The commutation relations of this subalgebra 
can be written simply by using the notation with the formal series. 

Lemma 2.3 S{u) obeys 

Qab 



[Saiu),Sb{v)] 



u — V 



+ 



Sa{u) - Sbiv) 



U + V 



(2.29) 



where Qab = ~ ^ab ^''^^ symmetry property that 



S{u) = S{-u)' . (2.30) 



Proof. The first of these is true by virtue of (j2.7p and the second is immediate from the definition 
([227]). ■ 

Note that these commutation relations can be obtained from the classical limit of the twisted 
Yangian introduced in [15]. More abstractly, the relations (j2.29p and (j2.30p can be regarded as 
defining an algebra, which can then be seen to be embedded in the half loop algebra according 
to (Il^TD. 

It is well-known that the centre of this subalgebra is generated by the odd coefficients of the 
series s{u) = tiS{u) (see for example [2T], section 4). But we have also 

Proposition 2.4 The quantities in the expansion of s'{u) = tr5(n)^ are mutually commuting, 
or equivalently 

[s'{u),s'{v)]=0 (2.31) 

for all values of u and v. Moreover, 

[5W,s'(n)] =0 . (2.32) 
The elements in S^^^ generate 5o{p, q) for rj = +1 and 5p{N) for rj = —1. 

Proof. The details of the proof are given in appendix [Bl ■ 
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3 Gaudin models 



3.1 The Inner-twisted Gaudin Magnets 

The quantum Gaudin magnet, introduced in [3j, is an integrable spin chain with long range 
interactions. The Gaudin Hamiltonians for the model with L sites are 

where Zi are complex numbers. (Recall that Pjk permutes the j*^ and k^^ spins.) This model 
is usually called the A^-type Gaudin model. It may be obtained from the more general class of 
integrable Hamiltonians 

Hk = j2 -^^^^^^ (3.2) 
j=i 

by specifying that the spin at each site j is in the fundamental representation of 

Now, given proposition 12.21 above, we can obtain new integrable models, as in the following 
proposition. These models describe spins placed at fixed positions in the plane, each of which 
interacts with the central point and with the other spins, not only directly, but also via their 
images under the rotation group of order n. 

Proposition 3.1 The model described by any one of the Hamiltonians 

=2:2: + E M 

is integrable. This model has gljy^ © gl^^ © • • • © Q^n^-i symmetry. 
Proof: From the definition (j2.18|) of B{u), one finds 

b'iu) = iTB{uf = jzYl ^\ 4"^ + E E (3.4) 
with H^^^ as given in the proposition. (The identity 



(3.5) 



(u — T 3zi){u — T ^Zp) T ^ Zl — T ^ Zp \U — T i Zl U — T ^ Zp 

for {j,l) 7^ ik,p) is helpful in showing this.) 

It then follows from proposition 12.21 that [Hj^\ Hp^^] = 0. Since (for n > 1) these operators 
Hp'^^ are independent we have found L commuting conserved quantities, completing the proof of 
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integrability of the Hamiltonian H^"'\ Next, from the relation [B^^\ trB{u)'^] = 0, also proved 
in proposition 12.21 we deduce that [B^^^H^"^^] = 0, which gives the Qljsf^ © glj^^ (£)■■■ (B 9^n„-i 
symmetry of the model. ■ 



Examples: 

• For n = 2, we obtain the Hamiltonian 

jj{2) _ / traPgfcFgj ^ tr aFakGaFajGa \ ^ traPqfcGalPafcG'a g-j 

of the BC-type Gaudin model studied in |19j . 

• If the sites carry the fundamental representation of gt^v, our Hamiltonian is 

Let us remark that in the Al case (n = 1) supplementary conserved quantities, called higher 
Gaudin Hamiltonians, can be found by computing for example trT(u)^ (see e.g. |22j)- The 
question of whether this is possible in the generalized cases (n 7^ 1) studied here remains open. 

3.2 The Outer-twisted Gaudin Magnets 

Using the algebraic result of proposition l2.4l we can also succeed in constructing integrable models 
based on outer automorphisms, as follows: 

Proposition 3.2 The model described by any one of the Hamiltonians 

^ Zk-Zj Zk+ Zj ) 2Zk 

is integrable. The model has 5o{p,q) symmetry (resp. sp{N) symmetry) for ry = +1 (resp. 

v = -i). 

Proof. The proof is similar to the one of proposition 12. 2i Using definition ()2.27p of S{u), we 
show that 



(3.9) 



9 



with given as in the proposition. Then, we deduce from proposition 12.41 that [Hj!,Hp] = 0, 
and since the operators Hp are independent for different p, this proves the integrabihty of H"^. 
The symmetry algebra is deduced from [S^^\ trS'(n)^] = proved in the proposition 12.41 ■ 

Every choice of representation Vi • • • V;, for the sites then yields a Gaudin-type model. (It 
is worth remarking that it is possible to choose different representations at different sites.) For 
example, in the fundamental representation of qIj^, the Hamiltonian is 




(3.10) 



We may interpret as a Gaudin model with boundary as in the BC type model (equation 
(j3.6p . and see also [19]). The Zk + Zj term in (j3.8p corresponds to the interaction between 
the k^^ spin represented in Vk and the j^^ 'reflected' spin transforming in the contragredient 
representation. This type of boundary is called soliton non-preserving and has been implemented 
in other integrable models |23t [2^ . The final term in (|3.8|) corresponds to the interaction 

between particles and the boundary. 

4 Calogero Models 

We turn now to the second class of integrable system of interest in this work, the Calogero 
models. We seek to construct dynamical models of multiple particles on a star graph, whose 
pairwise interactions are determined by a potential of the usual Calogero type, namely 
where q is the linear distance separating the particles in the plane of the star graph. We will first 
construct models of particles of unspecified statistics; subsequently, by specifying statistics and 
parity, we arrive at Calogero models for particles with internal spins. 

4.1 The Al case 

Let us first recall the Calogero model based on the root system Al [3], and in particular the use 
of Dunkl operators [27] in demonstrating its integrabihty [TB]. Consider a quantum mechanical 
system of L particles on the real line. Let qi be the position operator of the z*^ particle, and write 
the position-space wave function as 

ilj{qi,q2,...,qL)- (4.1) 
Let = be the operator which transposes the positions of particles i and j, 

..,qi,...,qj,...)=ij{...,qj,...,qi,...). (4.2) 
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Let us denote Sl the permutation group of L elements and (ij) the transposition of the elements 
i and j. Each element s € Sl can be written in terms of transpositions, namely s = (ij) . . . (kl). 
Then, we can define as the shorthand for the product i^ij . . . (even though the expression 
of s in terms of transpositions is not unique, is well-defined due to the commutation relations 
satisfied by l^ij). The sign of s, denoted is the number of these transpositions modulo 2. 

Let us define L operators di - the Dunkl operators - by [281 EH] 

1 d 

= Pi + A ^i,-, where -pi = -ih—. (4.3) 

Qi- Qj dqi 

It follows from the relations ^ijQj = Qi^ij that the Dunkl operators commute with one another, 

[di,dj] = (4.4) 

and consequently that the quantities 

JC^) = ^ (4.5) 

i=l 

form a commuting set also. The I^'^^ are algebraically independent for k = 1,2, ... L, and these 
give L commuting conserved quantities of the model with Hamiltonian 

i=l ^=l \ j^i J 

which is therefore, by construction, integrable. 

The next step is to consider particles with internal degrees of freedom, which we take to be 
in the fundamental representation of glj^. The wave funtion becomes 

Tp{qi,q2,.. ■,qL\vi,V2, . . ■ ,vl)- (4.7) 

where Vi G C''^. As we define operators which transpose the positions, we introduce operator 
Pij which transposes the spins 

Pijipiqi, ...,qL\---,Vi,...,Vj,...)= ipiqi, . . . ,qL\ . . . ,Vj, . . . ,Vi, . . .). (4.8) 

We define similarly to the matrix Ps = Pij . . . Pki for s = (ij) . . . (kl) acting on the spins. 

As explained before, to use the St Petersburg notation, we need supplementary spaces called 
auxiliary spaces (which are and, in this case, isomorphic to the quantum space) and denoted 
by the letters a, b,... The conserved quantities (j4.5p then emerge in a natural way from the matrix 



Ta(n)=X;-^, (4.9) 
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because (as one can see using tVaPai = 1) 

t{u) = tVaTa{u) 

Here (14. 9p is nothing but a modified version of tlie monodromy matrix (I2.10p . Tlie parameters 
zi are replaced by tlie Dunkl operators, and since the quantum spaces are chosen to be in the 
fundamental representation, P^^ = Eij (gi eji becomes the transposition operator Pai = Eij (gi Eji 
(for i = 1, . . . ,L). Now because the dj commute with each other and with all operations on the 
internal degrees of freedom, T{u) obeys the half loop algebra relations ()2.7p exactly as before. 

Suppose, finally, that the particles are in fact indistinguishable, which is often the case of real 
physical interest. One must then impose definite exchange statistics on the wavefunction: 

Pij^ijiP = ei;, (4.11) 

where e = +1 for bosons and e = — 1 for fermions. The projector onto such states is 

A=^el«l^,P, . (4.12) 

seSL 

The following relation 

(A-l)r(n)A = (4.13) 

demonstrated in [18J is crucial, because it implies that the modified generators T{u) = T{u)A 
preserve the condition Aip = ip, and obey the same algebraic relations as the original T{u). From 
T{u) we may define t{u) = iTaT{u)K = t{u)A, and hence = I^'^^A. Using [t{u),t{v)] = 0, 
one obtains that the I^^'^ are once more L commuting conserved quantities of the system with 
Hamiltonian 

where we are now able to replace which acts on particle positions, by P, which acts only on 
the internal degrees of freedom. Moreover, since t{u) commutes with T{u), the model has a half 
loop symmetry algebra. 

The subtlety in all this is that the Dunkl operators themselves do not obey any relation 
analogous to (j4.13p . There are thus essentially three steps in this procedure to construct an 
integrable Hamiltonian for a system of indistinguishable particles: 

1. Find commuting Dunkl operators, and hence T{u) 

2. Construct the appropriate projector A onto physical states, 

3. Prove the relation (A - l)r(n)A = 0. 




(4.10) 
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4.2 Dunkl Operators for the order n inner-twisted case 

We can now turn to applying these ideas to the model of interest in the present work. We consider 
a system of L particles living on n half-lines - "branches" - joined at a central node, as in figure 
[TJ The branches are given parametrically hy z = r^t, t > 0, A; G Z„, and we shall denote them by 

M+,rM+,...,r''-^M+. (4.15) 

As before, let qi be the position operator of ith particle. (Note that the spectrum of Qi is not 




Figure 1: A particle on the branch M^, and its images on the other branches 



real, but only for the superficial reason that we choose to regard the half-lines as subsets of the 
complex plane.) In addition to the ^ij, which exchange particle positions, we can define now 
new operators which move the particles between branches: 

^ii;i...,qi,...)=i;i...,Tqi,...). (4.16) 
It is useful to collect together the algebraic relations satisfied by the qi, cSj, and ^ij: 

CO). . CO) (M.. — ca., oi>. . ca., _ 1 cm. . — ca .. {AM\ 

= 1, (4.18) 

and 

^ij'^j — '^i^iji ^ij^j — Q.i^iji '^'^iqi — qi'^ii (4.19) 

with all the rest commuting. To construct an integrable model, the first task is to find a suitable 
generalization of the commuting Dunkl operators introduced above. 

Proposition 4.1 The Dunkl operators defined hy 

^-^ ^-^ qi — T'^Qi ' qi 




(4.20) 
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for arbitrary parameters A, /i^ G C, commute amongst themselves: 



[di,dj] = 0. 



(4.21) 



Proof. Consider first the terms at order A. We have 



(4.22) 



using the relations (|4.19|) and the definition pj = —ih-£p, which together imply JSjPj = pj^jT. 
The two terms of this type occurring in [di, dj] are 



1 



+ 



rt- G)k CO) . . O)- 



,,{T'p,-pj)^^^f (4.23) 



which cancel, after a change of the summation index in the second. The two terms containing 
[pj, l/{qi — qj)] cancel similarly. The terms occurring at order A'^ are of the form 



^ a}k cm. 6)-k ^ 



(4.24) 



_qi-'^(lg Qj-T'-qh 

These vanish trivially unless at least one of the indices i,g matches at least one of j,h. It is 
straightforward, though tedious, to check that the terms with exactly one index in common sum 



to zero, by using the relations (|4.19|) to bring every such term into e.g. the form ■^^-^^^^^JS.S 



J 1 

l-q q- 

and then summing the fractions directly. The terms in which both indices match give 



T ^qj - T^qi qj - r^qi ^qi - r ^qj qi - ^qj 



(4.25) 



and here the sum over (. may be re-written as 

1 

qi-T 



+ 



Qj - T^Qi Qj - r'^+^qi r ''qj - T^qi r '^qj - T^+^qi 



(4.26) 



which then vanishes by shifting the dummy index ^ in the second and fourth terms. The terms 
involving /i^ may be treated similarly. ■ 

These Dunkl operators have been introduced previously in ^30j as Dunkl operators associated 
to complex reflection groups. A proof of their commutativity is already given but is based on 
different computations. 



As in the Al case above, the quantities 



(4.27) 



i=l 
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are then mutually commuting, forming a hierarchy of Hamiltonians of an integrable system. Their 
detailed forms are rather complicated - for example, in the case of n = 3 branches with only 
L = 2 particles and jij- = 0, we find that the first three are 

/(I) = pi+p2 + Af ^1^-^+ ^r^^2^^i2 (4.28) 

\q1-rq2 qi-r^q2 ) 

+ \V2 (-^—^^2^1^ + ^12 
\q2-Tq1 q2-T'^qi J 

aQi + '^QiH + ^HQl + ql 



{qf-qiy 



7(3) = pI+pI (4.30) 

- 3Api (- ^(^12 + A) + L^(<^,2^-i^2 + A) + ?_^(^,2^,^-i + A)) 

V (qi - q2r iqi - ^q^Y {qi - r^q2r J 

- ^^P^ (1 — ^(=^12 + A) + - — L^(^,2^-i^, + A) + L_^(^i2^2^-i + A)) 

V [qi - q2r \12 - rqiY {q2 - T^qiV ) 

- 3V3iA2- ^— ^^,^2(^1 - ^2) 

\q\ - 92) (gi - Tg2)(gi - T^q2) 

4.3 Quasi-partity and peirticles with spin 

The next steps are to consider particles with spin and to choose a suitable projector onto physical 
states. We take the latter to be the product of two parts 

ApAq, (4.31) 

where 

kp=Y^ e\'W,Ps (4.32) 
is the projector onto states of definite exchange statistics, and 

= n ^ E (4.33) 

is a new projector which relates the wavefunction on different branches: = In the 

case of n = 2 this means relating the wavefunction at = a; with qi = —x, or, in other words, 
imposing a parity condition, so we shall call analogous conditions for arbitrary n "quasi-parity" 
conditions. The reason for requiring quasi-parity is that we would like to find a Hamiltonian that 
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does not involve S^ij or and, just as imposing definite statistics allowed 3^ to be replaced 
by P in (j4.14p . so here quasi-parity will allow =S to be replaced by G, the matrix defining the 
automorphism of 

Now Ta{u), defined in ()4.9p . does not respect quasiparity, but we have instead 
Lemma 4.2 

Ba{u) = r^GiTa{uT^)G-J (4.34) 

satisfies 

[5„(ii),AQ] =0. (4.35) 

Proof. By direct computation. ■ 
Thus, using also that T{u) obeys (|4.13p . we have 

(1 - KpKQ)Ba{u)kpKQ = 0. (4.36) 
and have therefore arrived at the following result 

Proposition 4.3 The modified generators B{u) = B{u)ApAq preserve the statistics and quasi- 
parity of the wave function and themselves satisfy the relations 12.20\) . The quantities I^^^ in the 
expansion of 

b{u) = traBa{u)ApAQ =:Y.^ (^.37) 

k=0 

are mutually commuting, and non-zero only when k = mod n. The model with this hierarchy 
of integrable Hamiltonians has symmetry g[j^[z]^ . 

Proof. Most of this follows from the construction above: it remains only to show that I^^^ 
vanishes for k ^ mod n. One sees this by writing 



e=i p=o \ fcgz„ 



^Note that in fact there is another natural class of quasi-parity condition: when n — 2, Tp{. . . , qt, . . . ) — 
Gitp{. . . , —qi, . . .) is obviously equivalent to tp{. . . ,qi, . . .) — Gi^{. . . , —qi, . . . ) = 0, but these two formulations 
suggest different generalisations to n > 2: we can on the one hand demand for all i that 

i>{. . . ,q^,. . .) = G,Tp{. . .,Tqr, ...) = ■•■ = G"~V(- • • ,T"~^'7i, ■ • ■) 

or, alternatively, for all i 

ip{...,qi,...) + rGiip(. . . , rqi, . . . ) H h r""^G"" V(- ■ • , T^'^Qi, . . . ) = 0. 

We use the first type of quasi-parity here. With such a condition in force one need only give the wavefunction 
on in order to completely specify the state, and in this sense the model is really on the half-line. The second 
type of quasi-parity is weaker - and so potentially interesting - but does not allow us to replace ^ by G in the 
Hamiltonian. 
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and noting that XlfeeZn unless p = mod n. ■ 

The vanishing of some of the generators is as expected: in the case n = 2, for example, the 
charge I^^^ , which is first order in momentum p, does not survive the introduction of a boundary. 
In the case of n = 3 with L = 2 particles, already mentioned in (j4.28p . the first non- vanishing 
charge is third order in momentum: 

7(3) = pI+pI (4.39) 

- 3Xpi (- ^-^(ePi2 + A) + ^{ePi2G^'G2 + A) + ^-^{ePuGiG^' + X)) 

- ^^P^ (7 ^(^^12 + A) + 7 r^iePuG^'Gi + A) + ^-^-^iePi2G2G^^ + A)) 

V [qi - q2Y (92 - rqiY (^2 - r^qiY J 

- 3-\/3zA^7 — — 5 — -GiG2{Gi — G2) ■ 

[qi - q2){qi - rq2){qi - T^q2) 

5 Conclusion 

In this paper, we gave the St Petersburg presentation of subalgebras of the gt^v half loop algebras 
associated to all finite order automorphisms of gt^y. This presentation allows us to obtain com- 
muting quantities used to prove integrability for new integrable models of Gaudin or Calogero 
type. The non-Abelian symmetry for each of these new models is also exhibited. 

We may expect that the usual methods to solve the Gaudin models and the Calogero models 
may be generalized to solve the models given in the propositions 13. H 13.21 or 14.31 introduced in 
this paper. Namely, for the Gaudin models, the Hamiltonians ()3.3p and (j3.8p may be diago- 
nalised by generalising the usual methods such as the separation of variables [31] or the algebraic 
Bethe ansatz [32|. For the Calogero models, the previous link established in [30] between the 
nonsymmetric Jack polynomials and the Dunkl operators (j4.20p may be useful to diagonalize the 
Hamiltonian given in proposition 14. 3[ 

Our discussion has dealt exclusively with quantum mechanical models. However, for each 
algebra introduced in the paper, there exists an associated Poisson bracket algebra, obtained 
simply by replacing the commutator on the left of the defining relations (j2.ip . (j2.3p . (j2.7p . (j2.20p 
and ()2.29p by a Poisson bracket. This allows us to treat certain classical mechanical problems. In 
such problems the entries of T{u), B{u) or S{u) are commuting functions on phase space, which 
simplifies many computations. For example, the results of propositions 12.21 and 12.41 are replaced 
by the stronger statements 

{bk{u),be{v)} = and {sk{u), se{v)} = (5.1) 

where bk{u) = tri?(n)'^ and Sfc(n) = trS'(ti)'^. These results strongly suggest that the classical 
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counterpart of models given by (j3.3p and (j3.8p are integrable in the sense of Liouville. (It remains 
to prove that the quantities are independent.) 

Finally, although the models of this paper were quantum-mechanical, the algebras are classi- 
cal, in the sense that they are not q'-deformed. A very interesting question is whether a similar 
construction of subalgebras from higher order automorphisms of gt^v is possible in the case of 
quantum groups. If so, then these subalgebras would be a g-deformation of those in this paper, 
and should also have associated to them integrable models on n half-lines. 

Acknowledgements: NC is grateful for the financial support of the TMR Network "EU- 
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HPRN-CT-2002-00325. CASY gratefully acknowledges the financial support of PPARC. 

A Proof of proposition 12.21 

It is convenient first to use 



G-'^PabG'^Baiu) = G-^GlB,{u)Pa, = GIb,{u){G,^ G^G'^ P^, = GIb,{u)G,^ G',^Pai,Gl 



[Ba{u),B,{v)] = y^{B,{v) + T^Ba{u)-T''Ba{r''v)-B,{T-^u)) " (A.l) 



keZn 

The goal here, and in the following, is to bring every term containing Pab into the form BaB^ G^^ PabG^. 
Next, we have 



and similar manipulations, to re-write the commutation relations (j2.20p as 



[Ba{u),Bb{vf] 



Bb{v) [Ba{u),Bb{v)] + [Ba{u), Bb{v)] Bb{v) 



(A.2) 




+ {Bb{v) + T^Baiu) - T^BaiT^v) - Bbir'^'u)) T^Ba{T%)} 



r^—k T3 rik 



u — t'^v 



{Bb{vf - Bb{v)Bb{r-^u) + r^^ Ba{u)Ba{T^v) - r^'^B^ir'^vf 



- T^BairK)Bb{T~''u) + T''Baiu)Bb{v) 

+ ABa{T'^v),Bb{T-''u)]-T''[Ba{u),Bb{v)]} 
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and then the brackets in the final hne may be evaluated by using (jA.ip once more, to give, after 
some manipulation of the summation indices. 



[Ba{u),B,{vf] = {Bb{v? -B,{v)B,{T-''u)+T^^Ba{u)Ba{T^v)-T^''Ba{T''vf 

keZn 

Ga '^PabGa 



'i3,(rS)i?fc(r-^) +r^i?„(n)i3,(z;) ° (A.3) 

J u — t'^v 

^T'^^Ba{T^v) - T'^^Ba{T^-^u) - T^Bk{T^-h 



U — T^-V 

j,kez 



+ r'^Bbiv) + T^+'Ba{u) - T^+''Ba{T^v)} 



{u — T^v){u — T^v) 



Therefore 



[Ba{uf,Bh{vf] = Ba{u)[Ba{u),Bb{vf] + [Ba{u),Bh{vf]Ba{u) (A.4) 

= Yl Ba{n){Bk{vf - B,{v)B,{t~'u) + T^''Ba{u)Ba{T^v) - T^^Ba{T^vf 



T^Ba{T^v)Bl,{T-^u)+T^Ba{u)B,{v)]'- 



a ^ab'-^a 
U — T^V 



+ Y [Bh{v? - B,{v)B,{t-^u) + T^^Ba{u)Ba{T^v) - T^^B^ir^vf 

ab 



r^—k p (^k 



U — T^V 



- T''Ba{T''v)Bb{T-''u) + T''Ba{u)B,{v)}T-''B,{T-''u)^ 

- Y Ba{u)[T^^Ba{T^v) - T^^ B^ir^-'^u) - T^B,{t^-^v) 
j,k€Z„ 

+ r'^B.iv) + T^+f^Baiu) - r^+'Bair^v)} . ^\ 'f'^ ' ^ . 

> yu — t'^v)[u — T-'v) 

- Y {r^^Bair^v) -T^^Ba{T^~''u) -r'^Bbir'^-^v) 
j,kez„ 

+ T^Bdv) + T^+^Baiu) - r^+'Bair^v)} ' Ba{u) 

J [U — T'^V){U — TJV) 

In the final line we could again add commutators to move all the BaS to the left of the -Bfe's, but 
in fact it is not necessary to do so in order to evaluate 

[h'{u), h'{v)] = t,ab [Ba{u)\B,{vf] . (A.5) 

Consider first those terms containing Pat- Because we avoided terms of the type Bj,BaBi,Pab in 
(|A.4p . these all reduce to single traces: for example 

tTabBa{u)Bb{vfG-''PabG'l = tl abG^M^) G' B,{vf Pab 
= tTabr-''Ba{T-''u)PabBa{vf = tTT~'' Ba{T-''u)Ba{vf , 
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where the final equahty is vahd because trbPab = 1- The remaining terms give products of traces, 
and, after some cancellation, one finds 



[b'{u),b'{v)] = ^Jll^tr 



fcez„ 



U — T'^V 



B{T-^uf,B{v) 



(A.6) 



+ { ^^r^B(u)B{T^v)G^-^ iT&-^ + tiB{T^v)B{u)G''-^ trG^" 



3r 



■i^B{u)G^-HiB{v)G^-'' - tiB{v)G^-^ trB{u)G^-^y 



(u — Tkv){u — T^v) 

The first term, cubic in B, reduces, given the identity tr[Af, A^] = tTah[Ma, Nb]Pab and ([Oil , to 
^ \^tTB{T^v)B{u)G^-HiG''-^ - tvB{u)B{T^v)G''-HiG^-'' (A.7) 



-T-^ tTBiu)^-^ tiB{v)G^-^ + iiB(v)a-^ t^B{u)G^-^\- 

) (U — T 

and, on collecting terms, one has 



[b'{u),b'{v)] = {r^'tr [B(T^■^;)G'=~^S( 



trG^" 



(u — t''v){u — T^v) ' 

(Ai 



+ 



tTB{v)G^-'',trB{u)G'''^ }- 



{u — T^v){u — T^v) 

The second commutator can be shown to vanish, and on evaluating the first one is left with 
[b'{u),b'{v)] = ^ {tr5(r-'n)G^-'trG'-^' - tr5(r-'u)G'-J trG'=-' (A.9) 

(r^' + r*^) trG^-'^ 



+ iiB{v)G^-HiG^'^ - iiB{v)G^-hTG^-^Y 



{u — T^v){u — t''v){u — T^v) 

Consider now the two terms containing B{v). After taking | the sum over the cyclic permutations 
of the dummy indices k, one finds that these reduce to 



iTB{v)G^-^iTG^-hTG 

j,k,l£Zn 



1-3. 



{u — T^v){u — t''v){u — T^v) 



and the coefficient of tr5(t;)G"°~nrG'' trG'' in this sum is (with a factor i when a = b) 



E 



+ 



j-h+l _ ^a-l 



(u — T"+'u)(n — r' "v){u — rv) (u — r"+'f )(u — r' °'v){u — t''v) 



(A.IO) 



(A.ll) 



which may be seen to vanish by using 



l/v l/v_ 



^ , and the same identity 



with (a ^ b). Similar arguments hold for the B(u) terms in ()A.9p . and we have, finally, that 

[b'{u),b'{v)] = 0. (A.12) 

It remains to show that = 0. This may be seen by expanding (|A.3|1 to leading order 

in 1/u and taking the trace in space b. The elements in B^^^ is the +l-eigenspace of a therefore 
they generate gl^^ qIn^^ © • • • © 9iN„-i- 
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B Proof of proposition 12.41 



Let us rewrite the relation (|2.29p as 



[Saiu), Sbiv)] = {Sa{u) + Sbiv) - 5fe(n) - Saiv)) 



ab 



+ 



U — V 



Saiu) - Sbiv), 



Q 



ab 



U + V 



iB.l] 



Note that, in contrast to the previous case computed in appendix|Al here the Q cannot be moved 
through the S. We are now in the position to compute the bracket [Saiu), Sbiv)'^] by again using 
(jB.ip to bring every term containing Pab on the right-hand side into the form SaSbPab- We find 



[Saiu),Sbiv)^] = ( Saiu) iSaiv) + Sbiv)) - Saiv)^ - iSaiv) + Sbiv)) Sbiu) + Sbiv)'' 



P 



ab 



U — V 
Qab 



+Saiu) Sbiv) - Sbiv) ^""^ Saiu) + iSaiu) - Sbiv)) Sbiv) , 

U + V U + V I U + V 



+ 



Saiv) - Sbiu) - Saiu) + Sbiv), 



VQa 



y2 _ y2 



(B.2) 



where we used also the property PabQab = vQab = QabPab- Now we can compute [S'a(u)^, Sbiv)'^~\ = 
Saiu) [Saiu), Sbiv)'^] + [Saiu), Sbiv)"^] S'a(n) and take the trace in spaces a and b. It is then 
straightforward to show that 

2 

-tv\S{u)\S{v]\ + 

U + V 



s'iu),s'iv)] = -^tr[Siu)\Siv)] + -^(tr[Siu)\Si-v)] - tr[5(-n)2, , (B.3) 

u — V U + V \ / 



where we have used the symmetry relation (12.30p and, for example, the following properties 

tlabSaiu)QabSbiu) = tTabSaiu)QabSaiu)^ = trS'('u)S'(n)'^ (B.4) 
i'CabQabSaiu)Sbiv)Saiu) = tl^abSbiu)^ Sbiv)Sbiu)'^ Qab = tr5'(M)^5'(v)5(ti)^ . (B.5) 



Next, using the property 

tr[5(a;)2,5(y)] = -tr[5(y),5(x)2] 



tr:[Saiy),Sbix)^]Pab 



(B.6) 



and relation we have that tr[5(x)2, = ^ tr[5(y), = ^ti[Saiy), Sbix)]Pab. 

Then using relation (g^H), we get tr[S'(x)2, = which implies that the R.H.S. of I^K^ 

vanishes and proves that = 0. Finally, expanding ()B.2p to first order in 1/u and 

taking the trace in space b yields [S^'^\trS'^] = 0. The elements in S'(°) is the +1 -eigenspace of T 
therefore they generate so(p, g) for rj = +1 and sp(A^) for r] = —1. This completes the proof of 
proposition 12. 4[ 
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